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Recently discovered fine structure of the morphotropic phase boundaries in highly piezoelectric 
mixture compounds PZT, PMN-PT, and PZN-PT demonstrates the importance of highly non-linear 
interactions in these systems. We show that an adequate Landau-type description of the ferroelec- 
tric phase transitions in these compounds is achieved by the use of a twelfth-order expansion of 
the Landau potential in terms of the phenomenological order parameter. Group-theoretical and 
catastrophe-theory methods are used in constructing the appropriate Landau potential. A com- 
plete phase diagram is calculated in phenomenological parameter space. The theory describes both 
PZT and PZN-PT types of phase diagrams, including the newly found monoclinic and orthorhom- 
bic phases. Anomalously large piezoelectric coefficients are predicted in the vicinity of the phase 
transition lines. 

PACS numbers: 77.80.Bh, 77.84.Dy, 81.30.Dz, 77.65.Bn 
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I. INTRODUCTION 



For many years, perovskite-type materials have been 
a subject of extensive research in both experimental and 
theoretical physics. On one hand, different representa- 
tives of the perovskite family exhibit a host of physical 
phenomena, such as piezoelectricity, ferroelectricity, and 
superconductivity; on the other hand, perovskite struc- 
ture is a relatively simple and, thus, attractive object for 
theoretical studies. 

Even though there is a long history of studies of 
perovskites, they still present new surprises. Recent 
X-ray and neutron diffraction studies on solid solu- 
tions PbZn_ x Ti a 3 (PZT), Pb(Mg 1/3 Nb 2/3 )i- :c Ti :c 03 
(PMN-PT), and Pb(Zn 1/3 Nb 2/3 ) 1 _ :!; Ti ;c 03(PZN- 

PTBfflffl have revealed new phases in the vicinity of 
the morphotropic phase boundary^ on the T — x phase 
diagram of the solutions. In a narrow Ti concentration 
range (x = 46 — 52%), the low temperature structure of 
PZT was found to be monoclinic (crystallographic 
symmetry Cm) with polarization vectot, P directed 
along [uuv],u < v pseudo-cubic direction^. A similar 
Ma structure has also been recently saen in PMN-PT 
below room temperature for x = 35%u. Another fine 
structure of the morphotropic phase boundary has been 
found in PZN-PT. For x = 9 — 11%, the low-temperature 
structure is orthorhombic (O, Amm2) with P || [101]u. 
Rhombohedral (R, Rim) unpoled crystal of PZN-PT 
(x = 8%) was also found to exhibit irreversible mono- 
clinic Mc (Pm, P || [Om;]) distortion when an electric 
field above a certain critical value is applied along the 
[001] pseudo-cubic direction! 

Early theoretical investigations of phase transitions in 
perovskites were concentrated on BaTi0 3 , which goes 
through a sequence of phases upon cooling: cubic(CL 
Pmim), tetragonal(T, P4mm), O, and R. Devonshire^ 
explains the behavior of BaTi0 3 within the framework of 



a phenomenological Landau-type expansion up to sixth 
order in terms of the ferroelectric order parameter - po- 
larization P. While successfully describing the phase di- 
agram of BaTi0 3 , the potential used in Ref. || lacks the 
high-order terms necessary to describe the low symme- 
try phases. Using a geometric argument based on the 
Curie principle, Zheludev and Shuvalova classified possi- 
ble positions of P with respect to the cubic unit cell. Due 
to the purely symmetric nature of this approach, it fails 
to distinguish between monoclinic phases Ma and Mb 
(P || [uuv],u > v), because these phases have the same 
crystallographic symmetry Cm. The group-theoretical 
relationship between the geometric method and Landau 
approach was established in Refs. [ToUll| . 

Using this approach, Gufan and SakhnenkcEl found 
that on a two-dimensional (e.g. T — x) phase diagram 
of perovskites there can be a point (Tq,Xo) where five 
phases C, R, O, Mc, and T co-exist. They calculated 
the phase diagram in the vicinity of this five-phase point. 
Only small solutions of the equations of state that are 
close to the five-phase point were considered, justifying 
expansion in powers of small parameters (T— To), (x—xq). 
However, the results of this work do not apply to the 
recently discovered phases of the mixture compounds of 
Pb-based complex oxides since, in this case, the R-O-T 
and R— M^— T triple points are separate from the C-R-T 
triple point. Therefore, theoretical consideration cannot 
be limited to small solutions of the equations of state, 
especially for the lowest symmetry phases. 

Ab miiioEH, as well as phenomenological0, calcula- 
tions, have been used to account for the presence of mon- 
oclinic phases on the T — x phase diagrams of ferroelec- 
tric perovskites. Vanderbilt and CohcnEJ calculate the 
phase diagram in the space of phenomenological parame- 
ters within the framework of Landau-Devonshire theory. 
Their model is based on the eighth-order expansion of 
the Landau potential in terms of the polarization orien- 
tation P/|P|. Although monoclinic phases appear in the 
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phase diagram of the model, it does not incorporate cubic 
and triclinic (Tri, PI) phases. 

Having considered a number of successively more com- 
plicated models, a natural question arises: What is the 
most general phenomenological model of the phase dia- 
gram of a cubic system induced by a ferroelectric order 
parameter? This question can be answered by use of the 
concept of integrity rational basis of invariants (IRBI), 
introduced in Appendix A. IRBI can be thought of as a 
basis in the space of polynomials (formed from the order 
parameter components), which are invariant under the 
transformations of the symmetry group of the system. 
By a group-theoretical argument, if the IRBI contains 
polynomials of maximal order n, then at least a 2n-th- 
order phenomenological model is necessary to describe 
all the possible phases induced by the order parameterlia. 
This statement is true for the irreducible representations 
of groups generated by reflections, including the Pm3m 
symmetry of the perovskite structure. We will show that, 
in the case of perovskites, the Landau potential has to be 
expanded up to the twelfth-order terms to describe the 
phase diagram induced by the ferroelectric order param- 
eter. 

The results of the analysis of a simple twelfth-order 
model are presented in this paper. In Sec. II, we present 
the solutions resulting in a phase diagram containing all 



the phases allowed by the symmetry of the ferroelectric 
order parameter. In Sec. M , we fit the T — x phase di- 
agrams of PZT and PZN-PT. In Sec. |v|, we discuss the 
piezoelectric properties of the compounds in the vicin- 
ity of the newly found phase boundaries. Appendices 
are intended to give a group-theoretical and catastrophe- 
theory background for some general statements we make 
in the text. 



II. TWELFTH-ORDER LANDAU-TYPE MODEL 

Three algebraically independent Pm3m-invariant 
polynomials can be formed (see Appendix A) from the 
polarization vector components (P x , Py, Pz)- 

■h = PlPl + P^Pl + PlPh (1) 
J 3 = PlPlPl. 

The ferroelectric part of the Landau potential can 
then be expressed in terms of algebraic combinations of 
J i, J2, and J3. Since polynomials of up to sixth order 
are present in the IRBI ([[]) , the Landau potential has to 
be expanded up to the twelfth order terms: 



a\ Ji + a 2 Ji + b\ J 2 - 2nd and 4th order terms 

+03 Jf + di2 Ji J2 + c\ J3 ~ 6th order terms 

+a^J^ + (iii2^i J-2 + &2^f + disJiJa - 8th order terms (2) 

+a 5 Jf + dnx2 Jf Ji + di22</i Jf + dn 3 Jf J 3 + d 23 J 2 J 3 - 10th order terms 

+a 6 Jf + dmi2 Jf J2 + dii22 J\ Jf + hJ 2 + dim Ji J3 + ^123^1 J2 J3 + c 2 Jf - 12th order terms 



A complete investigation of extrema of the function 
F(P) in the multi-dimensional parameter space is a 
rather tedious exercise. However, main features of the 
phase diagram can be obtained from simplified models 
based on potential (Q) with some terms omitted. The 
truncated potential should satisfy at least two require- 
ments: i) it has to be bounded from below, and ii) 
it should, be structurally stable in catastrophe theory 
senseEj'tj. The latter requirement provides that small 
perturbations, which can arise from the omitted in (j^) 
terms, do not drastically change the results obtained in a 
simplified model. In Appendix [b], we briefly describe how 
compliance with the second requirement can be verified. 

A twelfth-order model 

Fvait = a\J\ + biJ 2 + ci J 3 + a 2 Ji + h Jf + c 2 Jf (3) 

meets both requirements when a2,b 2 , and C2 are positive 
quantities, while <Xi,&i, and c\ are parameters driving 
phase transitions, and can be of any sign. In spite of its 
simplicity, model (0) is in full agreement with the results 



of the group-theoretical analysis of the ferroelectric phase 
transitions in perovskites. It gives an account of all the 
phases C, T, O, R, Ma, Mb, Mc, and Tri, and it correctly 
describes the phase boundaries. For example, although 
the symmetry groups of the phases M^(Mb) and R obey 
a group-subgroup relation Cm C i?3m, phase transitions 
R-Ma and R-Mb cannot be of second ordeiOEJ. 

A complete phase diagram of the Landau potential (|^) 
can be constructed in the space of phenomenological pa- 
rameters (<2i,6i,ci). All the characteristic features can 
be seen from the two-dimensional cross-sections in the 
aioi-plane, as shown in Fig. [I]. 

Below we show how the phase diagram of Fig. [I] was 
obtained, using the example of phase R. Minimizing the 
potential (||) with respect to P x , P y , P z , and then impos- 
ing the condition P x = P y = P z = P s /y/3, we obtain the 
value of spontaneous polarization P s as a solution of the 
equation of state 
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(b) (c) 



FIG. 1: (a), (b) Phase diagram in ai&i-plane for the potential Fmu for c\ < 0(a) and ci > 0(b). Solid lines - first-order phase 
boundaries; dashed lines - second-order phase boundaries; dotted lines - stability boundaries of phases, (c) Diagram of possible 
second-order phase transitions between the phases. 
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The parametric equations for the boundaries of the 
phase stability domain are obtained by replacing the in- 
equality in either (|) or (|^) with an equality and solving 
it together with (W). The line resulting from (|) and (§) 



has a cusp at c\ < shown in the top inset of Fig. [j](a) 
This feature is defined by the parametric equation 



ci + 126 2 P 



cusp 



20C 2 a 

-P e = 

2J CUSp 



(8) 



This equation has no real solutions for c\ > 0, as the cusp 
merges with the a\ axis. Another interesting feature of 
the phase diagram is a first-order transition line below the 
b\ axis, shown in the top inset of Fig. [l](a) . It corresponds 
to an isostructural transition between two phases with 
the same structure and symmetry R, with the value of P s 
being the only difference between themO. Let P' s and P'J 
be the two values obtained from the equation of state (Q). 
Then the transition line is defined by Fr(P' s ) = Fr(P"), 
where Fr(P s ) is the potential of phase R: 



ai(P, 
h(P s , 



9 s 9 s 243 s 
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(7) 
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(9) 
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All the equations relating to isostructural phases have 
an additional permutation symmetry Pj 2 <-> P'J 2 . In- 
troducing new variables U = (Pj + P'J )/3 and V = 



P' s PJ /9, which are invariant under the permutation, 
we obtain 



6ai - Cl (U 2 + 2V) - 18b 2 U 3 - 2c 2 U 3 (2U 2 - W) = 0, 

ax - 4(3a 2 + b x )U - c x {2U 2 + V) - l2b 2 U(2U 2 - V) - 2c 2 U 3 (2U 2 - 5V) = 0, (10) 
2(3a 2 + h) + c x U + 12b 2 (U 2 - V) + 2c 2 (U 4 - 3U 2 V + V 2 ) = 0. 



V can be eliminated from ([To]) and resulting equations 
can be solved for a x and b\. One of the solutions para- 
metrically defines an isostructural R-R transition line by 



where i is a real parameter. The coordinates of the 
four-phase R-Ms-M^-Tri critical point can be obtained 
from (nT 



ai{U) = 

h(U) = 



18c 2 U 

„2 
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F36 2 f/ 2 - 
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(11) 

A similar isostructural transition line appears in phase 
Ma and is shown enlarged in the bottom inset of Fig. [I]. 

The first-order R-C transition occurs when the two 
phases have equal potentials: Fr(P s ) = 0. This condi- 
tion together with equation of state (0) gives 



ai(P s ) 



P 4 



-**I* 
9 

bl (P s ) = -3a 2 -^P 2 - 
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(12) 



We would like to point out some other important fea- 
tures of the phase diagram of Fig. 0, which can be rep- 
resented in a simple analytical form. For c\ > 0, phase 
Mc lies between T-Mp and O-Mc second-order phase 
boundaries defined respectively by 



bi = 0, and b\ 



a 2 b 2 



(13) 



The second-order O-M^ transition line is parabolic in 
the (ai,6i) parameter space 



(26 2 ai + ci&i) 2 - (4 + 8a 2 b 2 )(ciai - 4a 2 6i) = 0. (14) 

The lowest symmetry Tri phase appears at c\ < be- 
tween phases and Mb in a slice restricted by a para- 
metrically defined line 



ai(t) = -4a 2 t 2 + 
h{t) = -2b 2 t 4 + 



cia 2 

c 2 t 4 ' 
2gb 2 

c 2 t 2 > 



(15) 



a-i = 3a 2 



4 Cl 
C2 



1/3 



1/3 



(16) 



The first-order transition lines M B -R, R-M^, Ma-T, 
and R-O, and isostructural line were obtained 

numerically as the method described in Eqs. (^)-(p^|) in 
these cases leads to high order equations for a\ and b\. 

All possible second-order phase transitions are shown 
in Fig. 0(c). In agreement with the general analysis^, 
only phases R and T are accessible from phase C via a 
second-order transition. The M^-Tri second order tran- 
sition occurs in the c\ = plane and is not shown in 

Fig. g. 

In Sec. |n] we will give special attention to some other 
features of the phase diagram shown in the lower inset of 
Fig. |l](a). They are related to the structure of the phase 
diagrams, recently observed in PZTuu and PMN-PTO. 



III. PHASE DIAGRAMS OF PZT, PMN-PT 
AND PZN-PT 



The seemingly simple twelfth-order model potential 
Finit from Eq. (||) has led, as can be seen from Fig. [I], 
to a complex phase diagram describing all eight phases 
allowed by the symmetry of the problem. All pos- 
sible second-order phase transitions among the proper 
ferroelectric phases in perovskites are established by 
Fig. ^(c) , so any approximations made to describe transi- 
tions driven by the ferroelectric order parameter should 
not add or forbid any second-order transition lines ob- 
tained in the model (||). On the other hand, first-order 
phase transitions between any two of the eight phases 
can be expected, since, in general, there are no argu- 
ments that can forbid these transitions. For example, a 
direct O-Tj-transition has been experimentally observed 
in PZN-PT13 and earlier in BaTi0 3 . 

In Fig. ^| we plot the experimental phase diagrams of 
PZT and PZN-PT together with our calculations based 
on the model described. In Sec. II we have shown that 
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FIG. 2: (a), (c) Phase diagrams of PZT and PZN-PT from Refs. ^J§[18|. (b), (d) Phase diagrams calculated based on the 
model potentials Fpzt and Fpzn, as described in the text. Critical points are shown as triangles, their coordinates are marked 
on axes. 



the triple point of phases R, T, and Ma can appear as 
an intersection of the first-order phase transition lines 
(see Fig. 0(a)). A second-order phase transition M^-T 
is possible below the critical point S. This is exactly the 
topological structure of the MtT phase transition line 
that has been observed in PZTtil. To obtain reasonable 
agreement with the experiment, we assumed that the co- 
efficients ai and bi in (^) are linear functions of temper- 
ature and composition, while the rest of the parameters 
are constant: 



a\ = a T (T - T ) + a x (x - x ), , 
h = (3 T (T-T )+l3 x {x-x ). [L '> 

Keeping in mind the global topology of the phase di- 
agram shown in Fig. |l|, we can identify the features re- 
lated to the piezoelectric compounds. When applying 
the developed model to the particular systems, we add 
some symmetry allowed terms to the expansion (||) . This 
enables us to use simple temperature and concentration 
dependencies of the phenomenological parameters while 
achieving a good fit to the experimental phase diagrams. 
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TABLE I: Parameters of the Landau potentials used for qual- 
itative fit of phase diagrams for PZT and PZN-sPT. 



T- ' iA r a tyi of or 
JT aL dlllc LcI 


vdiLic iui r PZT 


V dlLlc IUI -T P Z N 


Q T 


1 
1 


1 
1 




-U.U J 


-u.o 


Pt 


u.uo 


n « 


Ac 


i 


t 


Cl 


-25 


-40 


a 2 


9 


9 


hi 


0.8 


0.7 


C2 


3.6 





dl2 


-0.3 


8 


di3 





13 



To model a PZT-type phase diagram (also related to 
PMN-PT), we use an expansion for the Landau potential 

FpZT = F init + d\ 2 J\J 2 . (18) 

Upon this modification, the M^-T second-order transi- 
tion line remains straight but becomes tilted with respect 
to the ai axis: 



Since the phcnomcnological phase diagram can always 
be isomorphically transformed to fit the experimental 
data, we do not try to obtain a perfect quantitative fit. 
Rather, we calculate the coefficients in arbitrary units to 
reproduce the qualitative features of the phase diagrams. 
The calculated values are shown in Table |. 

We used a relatively large absolute value for c\ to ob- 
tain comparable lengths of the first-order M^-T transi- 
tion line and R-T phase boundary, modelling the exper- 
imental data for PZT. This results in a vanishingly small 
region of coexistence of the two isostructural phases 
(see the bottom inset in Fig. 0(a)). The critical point 
S then almost merges with point L, the latter being a 
tangency point of the T-M^ second-order transition line 
continuation and the stability boundary of phase M^. 

For the expansion (Uq), point L has coordinates: 



We find that the part of the global phase diagram, 
relevant to PZN-PT, is almost unaffected by the variation 
of c 2 . Therefore, we assume c 2 = 0. The other coefficients 
are found in Table I. 



IV. PIEZOELECTRIC PROPERTIES 

Elastic energy and electrostrictive coupling have to be 
introduced into Landau potential to obtain the compo- 
nents of the piezoelectric tensor dik- The piezoelectric 
constants are inversely proportional to combinations of 
the second derivatives of Fpzr and Fpzn with respect 
to P x ,P y , P z cB. Some of these combinations vanish at the 
stability boundaries of the phases, resulting in anomalous 
values of the piezoelectric constants. 

By analyzing the phase stability conditions obtained 
in Sec. II, we find the components dik which can be large 
for each phase of the phase diagrams Fig. |[ 

In the R phase, all four independent piezoelectric con- 
stants are expected to be large at all the phase transition 
lines with phases T, M^, and O. 

The tetragonal piezoelectric dis (in this section we use 
a pseudo-cubic orthogonal basis in the dik notation) com- 
ponent of the piezoelectric tensor can be large at the T- 
R, T-M^ and T-0 phase boundaries, while (J33 and (J31 
are not expected to have anomalies in temperature and 
composition dependence at these lines. 

A complex behavior of the piezoelectric constants is 
expected in phase Ma at the boundary with phase T. In 
the second-order transition part of this phase boundary, 
no anomalies are expected. This is because no stabil- 
ity condition is violated there, even though it is a bifur- 
cation line where the monoclinic solution of the equa- 
tions of state becomes imaginary. However, this situa- 
tion changes drastically when M^-T is a first-order phase 
transition. In this part of the M^-T boundary and at the 
M^-R boundary, all ten monoclinic piezoelectric coeffi- 
cients can be large. 

Another verifiable prediction of our theory is the con- 
nection between the proximity of phase to the tri- 
clinic distortion and anomalies in the monoclinic piezo- 
electric coefficients dn, di2, di4, and dis, at low temper- 
atures. 

Finally, d33 and d3i can be large in phase O at the T 
boundary, and d24 - at the R boundary. 



ax 



4a 2 &2 - d{ 2 



^a 2 b 2 



d\ 2 



(20) 



More extensive modifications of Fi n a have been used 
to accommodate the features of the experimental phase 
diagram of PZN-PT, namely a direct O-T transition and, 
more importantly, the separation of the triple points R- 
O-T and C-R-T: 



F. 



PZN 



^12 Ji J2 + disJiJs 



(21) 



V. SUMMARY AND DISCUSSIONS 

We have shown that the increasingly high-order phe- 
nomenological models, used to describe the highly non- 
linear piezoelectric systems PZT, PMN-PT, and PZN- 
PT, are approximations of a complete twelfth-order 
model (0). Such a complex model is difficult to ana- 
lyze, but even the significantly simplified model poten- 
tials Fpzr and Fpzn describe the phase diagram of the 
compounds reasonably well, as seen in Fig. 0. We have 
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allowed for some discrepancy between the calculated and 
the experimental phase diagrams, for example, in the lo- 
cation of R-Ma and R-0 phase boundaries, in order to 
keep the dependence of the phcnomcnological parameters 
on temperature and concentration as simple as possible. 
On the other hand, it is clear, that since our phcnomcno- 
logical models correctly describe the qualitative features 
of the phase diagram, a complete quantitative agreement 
can be achieved by a careful fitting of the phenomenolog- 
ical parameters to the experimental data. 

The important features that our simple models cor- 
rectly reproduce are the separation of the triple points 
R-O-T and R-M^-T from the C-R-T triple point, and 
the switching of the M^-T phase transition between the 
first- and the second-order types. 

We also predict an anomalous behavior of certain 
piezoelectric constants along the phase transition lines. 
These predictions elucidate the nature of the ultra-high 
piezoelectricity and make our theory easily verifiable. To 
our knowledge, systematic measurements of evolution of 
the piezoelectric constants in the vicinity of the newly 
found phase boundaries in PZT, PMN-PT, and PZN-PT 
have not been done yet. 

Finally, group-theoretical and catastrophe-theory 
methods of phase transition theory described in this pa- 
per stand on their own. We described a general approach 
that puts phenomenological analysis on firm footing, and 
minimizes the necessary amount of numerical calcula- 
tions. 



where each Ci is an invariant subgroup of A+i and E - 
identity group. 

The problem of constructing IRBI of L is reduced to 
constructing invariants of the factor groups 

Ai = A/A-i (A2) 

at every step of ( |Al|) . It turns out that for any £-group, 
Ai acts on the invariants obtained in the previous — 
th step as cyclic group of order 2 or 3, for which the 
invariants can be easily constructed. Removing the al- 
gebraic dependencies among thus obtained polynomials 
at each step, at N-th. step one obtains the IRBI for the 
given L- group. 

If m is the number of invariants, constituting the IRBI 
for a space group representation then, in general, m > p, 
and there exist m — p algebraic relations of higher than 
first order among the m invariant polynomials. However, 
for irreducible representations of groups generated by re- 
flections m = p. This is the situation for the ferroelectric 
order parameter in Pm3m group, for which the step-by- 
step calculation of invariants results in the IRBI (Q) . 

The method described here allows for a simple gener- 
alization to groups with continuous subgroups, such as a 
gauge transformation groupOO. 

APPENDIX B: STRUCTURAL STABILITY OF 
PHENOMENOLOGICAL POTENTIALS 
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APPENDIX A: CALCULATION OF IRBI FOR 
REPRESENTATIONS OF SPACE GROUPS 

Let C be the group of all different matrices of a p- 
dimensional representation of a space group. The C- 
group can be thought of as a point group in the space 
of the order parameter components i = 1, . . . ,p. The 
Landau potential is invariant under transformations of 
this group. 

Integrity rational basis of invariants (IRBI) of the C- 
group is a minimal set of homogeneous /^-invariant poly- 
nomials depending on rji, such that any other invariant 
homogeneous polynomial is an algebraic combination of 
the elements of this set. O-ae of the methods to find IRBI 
is described briefly belowtil. 

It can be proved that for every group C the following 
sequence of subgroups can be constructed: 



In this appendix we review the method based on catas- 
trophe theory, which is used to analyze the stability of 
the potential functions. We consider the simplest case 
of an irreducible representation of a group generated by 
reflections, for which m = p (see Appendix 0) . 

To determine the structural stability of some potential 
function F(J\, . . . , J m ), we need to introduce the alge- 
braic combinations: 



Ui(Jx, . . . , J m ) 



fe=i K 



1, 



£e£„c£iC...C/; 



N 



c. 



(AI) 



( B1 ) 

where VJi is the gradient of invariant Ji in the order 
parameter space. Scalar products that appear in (Bl) 
can always be expressed in terms of invariant polynomi- 
als Ji, . . . , J m . A term can be omitted in the potential 
function without violating the type of the extremal be- 
havior of this function if its coefficient is small and if it 
can be expressed as 



^V i (J 1 ,...,J m )U i (J 1 ,...,J m ) + h. o. t. (B2) 

i=i 

Here Vi(J\, .. . , J m ) are some polynomials and 'h. o. t.' 
stands for 'higher order terms'. 

Following this algorithm, we obtain for the integrity 
basis (0) and potential ([|) 
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It is now clear that every term i n (B ), additional to (||), 
can be represented in the form (|B2| ). However, there is 
(VJi) = 4Ji, (V Ji, V J2) = 8J2, another restriction that arises from the requirement for a 

(VJ2) 2 = 4Ji J 2 + 12J 3 , (V J2, V J3) = 8J1 J3, (B3) model to produce all the phases allowed by the symmetry 
(VJ3) 2 = 4J2^3, (VJi, VJ3) = 12J 3 . 01 the order parameter. As can be seen from Eq. jl^), 

triclinic phase cannot appear in the phase diagram of the 
id truncated below the twelfth-order Landau potential ex- 

pansion (||). Consequently, the potential (g) is the sim- 
XJ\ = 4ai J\ + h. o. t. plest model that meets all the requirements mentioned 

U 2 = 8ai J 2 + h. o. t. above. 
U 3 = 12<2i J 3 + h. o. t.. (B4) 
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